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Confinement has been suggested as a tool to tune the self-assembly properties of nanoparticles, 


surfactants, polymers and colloids. In this way, we explore the phase diagram of Janus nanoparticles 


confined between two parallel walls using molecular dynamics simulations. A nanoparticle was modeled 


as a dimer made by one monomer that interacts via a standard Lennard Jones potential and another 
monomer that is modeled using a two-length scale shoulder potential. This specific design of the 


nanoparticle exhibits distinct self-assembled structures and a water-like diffusion anomaly in the bulk. 
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Our results indicate that besides the aggregates observed in bulk, new structures are observed under 
confinement. Also, the dynamic and thermodynamic behavior of the fluid phase is affected. The systems 
show a reentrant fluid phase and density anomaly. None of these two features were observed in bulk. 


Our results show that geometrical confinement leads to new structural, thermodynamical and dynamical 
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Introduction 


Colloids and molecules with anisotropic shapes and interactions 
play a significant role in condensed matter physics, especially in 
the design of self-assembled structures.”” In particular, Janus 
colloids are characterized as particles composed by at least two 
physically or chemically distinctive surfaces. They can have 
different shapes, as rods, spheres and dumbbells. These systems 
have a wide range of applications including medicine, self-driven 
molecules, catalysis, photonic crystals, stable emulsions, bio- 
molecules and self-healing materials.” 

Among the distinct Janus particle shapes, Janus dumbbells 
are colloids formed by two spheres, each one with distinct charac- 
teristics linked together with a separation that varies from an 
almost total overlap to one or two monomer diameters. Due 
to the resemblance between Janus particles and competing 
interaction systems, Janus dumbbells behave as surfactants in 
water-based emulsions due their amphiphilic properties.’® ** 
Self-assembled lamellar or micellar phases were observed in these 
systems due the competition between attractive and repulsive 
forces. 

Recent studies reported the production of silver-silicon 
(Ag-Si)" and silica-polystyrene (SiO,-PS)”* Janus dimers. Silica 
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behavior for this Janus nanoparticle. 


and silicon are classified as anomalous fluids, and therefore 
have a set of properties that diverge from the observed for 
regular materials. For most part of the fluids, the diffusion 
coefficient decreases when the pressure (or density) increases. 
However, materials such as water,”° silicon” and silica’? show 
diffusion anomaly, characterized by a maximum in the diffusion 
coefficient at constant temperature. Besides diffusion (or dynamical) 
anomaly, water, silicon, silica and others fluids, the so-called 
anomalous fluids, also have other classes of anomalies, such as 
structure and thermodynamic anomalies. In particular, the 
density anomaly is characterized by the increase of density 
with temperature at a fixed pressure. 

Distinct computational models were proposed to study the 
anomalous behavior of fluids. Among these models, an effec- 
tive two length scale (TLS) core-softened shoulder potential 
is an interesting tool to investigate systems with water-like 
anomalies. In particular, the model proposed by Oliveira 
et al.”*° reproduces qualitatively the diffusion, structural and 
thermodynamic anomalies, and was broadly used to study 
anomalous systems. This effective approach to describe anom- 
alous fluids was used to study monomeric and dimeric systems 
of anomalous particles.”° *? The TLS potential was used in our 
previous works to study the behavior of Janus dumbbells 
composed of one anomalous and one non-anomalous monomer 
in bulk.***° Despite the presence of the non-anomalous monomer, 
the diffusion anomaly was preserved. 

Confinement was proposed as an approach to tune the self- 
assembled morphologies. By controlling the confinement 
intensity it is possible to create micelles with distinct shapes. 
Computational and experimental studies have already explored 
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the confinement effects of the self-assembly of polyhedral nano- 
particles,*° patchy spherical colloids,*” asymmetric and sym- 
metric dumbbells,***° surfactants and polymers.*°** Moreover, 
confinement affects the diffusivity of spherical Janus swimmers." 
In fact, confinement strongly affects the behavior of fluids in 
general. For the case of anomalous fluids,“ new anomalies can 
be observed due confinement,” and even a superdiffusive regime 
can be induced.“ The anomalous region in the pressure x 
temperature (PT) phase diagram is usually shifted due to confine- 
ment. This shift can be to higher or lower temperatures, based on 
the nature of the fluid-wall interaction.”° 

Therefore, the question that arises is how the confinement 
will affect not only the self-assembled structures, but also the 
dynamical and thermodynamical behavior of the anomalous/ 
non-anomalous Janus dumbbell system. In this way, we perform 
intensive Molecular Dynamics (MD) simulations of Janus nano- 
particles composed of anomalous and non-anomalous mono- 
mers confined between two flat plates. In addition to the myriad 
of structures tuned by confinement, we show how the confine- 
ment and the TLS potential lead the system to have not only 
diffusion anomaly, but also density anomaly, not observed in 
the bulk. 

The paper is organized as follows: first we introduce the 
model and describe the methods and simulation details; next 
the results and discussion are given; and then we present our 
conclusions. 


The model and the simulation details 


In this paper all physical quantities are computed in the 
standard Lennard Jones (LJ) units,*” 


* 


= 


, p’=po, and =i (1) 


; 
o mo? 


for distance, density of particles and time, respectively, and 


p =— and T* =— (2) 


for pressure and temperature, respectively, where o, € and m are 
the distance, energy and mass parameters, respectively. Since 
all physical quantities are defined in reduced LJ units, the * is 
omitted, in order to simplify the discussion. 

The systems have N = 1000 dimeric particles, totalling 
N = 2000 particles, confined between two smooth and parallel 
plates. The Janus dumbbell particles were modeled using two 
spherical core-softened particles, each one with mass m and 
effective diameter c, linked rigidly at a distance 4 = 0.8. The 
dimers are formed by type A and type B monomers. 

The particles of type A present anomalous behavior and 
their interaction is given by a two length scale potential, 
potential AA, defined as*”*° 


une) 7 ‘|(2) Š (=) ] + uo exp] 4 (2 = o » (3) 


where rj = |7; — 7;| is the distance between two A particles i and j. 
The first term of the potential is a standard 12-6 LJ potential*” 
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and the second one is a Gaussian shoulder centered at 79, with 
depth uo and width co. The parameters used in this work are 
Up = 5.0, Co = 1.0 and r/o = 0.7. Both systems, monomeric and 
dimeric, modeled by this potential, present density, diffusion and 
thermodynamic anomalies, like those observed in water, silica and 
other anomalous fluids.” +4849 Once the monomeric system of 
A particles shows water-like anomalies, we will refer to this specific 
monomer as an anomalous monomer or an anomalous particle. 

The interaction between particles of type B, the potential BB, 
is given by a standard 12-6 LJ potential, like the first term of 
eqn (3), cut and shifted at the cutoff radius re, 


Urs (ry) — Urs(re), fy < fe, 
UM (14) = (a) 


0, Fij > Fo. 


The BB potential has a cutoff radius of r, = 2.5. Meanwhile, 
the interaction for A-B particles is given by the Weeks-Chandler- 
Andersen (WCA) potential, defined using eqn (4) with r, = 2" as 
the cutoff. The interactions between dimers and walls are given by 
the projection of the WCA potential in the z-direction. The 
potentials are illustrated in Fig. 1. 

The simulations were performed in the canonical ensemble 
using the ESPResSo package.°”°* The number density, defined 
as p = N/V, where V = L” x L, is the volume of the simulation 
box, was varied from p = 0.05 to p = 0.50. In all simulations, 
L; = 4.0 and L was obtained from L = [N/(pL,)|"”. This separation 
Lz was chosen to simulate the case of a strong confinement, 
where the system is quasi-2D. Standard periodic boundary 
conditions are applied in the x and y-directions. The tempera- 
ture was simulated in the interval between T= 0.05 and T= 0.60. 
The system temperature was fixed using the Langevin thermo- 
stat with y = 1.0, and the equations of motion for the fluid 
particles were integrated using the velocity Verlet algorithm, 
with a time step t = 0.01. We performed 1 x 10° steps 
to equilibrate the system. These steps are then followed by 
5 x 10° steps for the resulting production stage. To ensure that 
the system was equilibrated, the pressure, kinetic energy and 
potential energy were analyzed as a function of time, and 
several snapshots were taken at distinct simulation times. 


Sr e e oe 
F *—- Potential AA 
4r ih Potential AB - 
--—- Potential BB 


Fig. 1 Interaction potential between particles of type A (dot-dashed blue 
line), between particle of type A and B (solid magenta line) and between 
particles of type B (dashed red line). The interaction between dimers and 
confining walls is given by the projection of the potential AB in the 
z-direction. Inset: Janus dumbbells formed by A-B monomers. 
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Since confined systems can be sensitive to the number of 
particles in the simulation, for some points we carried out 
simulations with 5000 and 10000 particles, and essentially the 
same results were obtained. Moreover, we run some points with 
a production time of 1 x 10° to test if the system was well 
equilibrated, and the same results were obtained. 

The system dynamics was analyzed using the lateral mean 
square displacement (LMSD) as a function of time, given by 


(icm) = Ficm(to)|’) = (Aficm(t)), (5) 


where Fiem(to) z (Xem(to)” a Yemn(to)’) and Ficm(t) E (elt) $ Yem(t)”) 
denote the parallel coordinate of the nanoparticle center of 
mass (cm) at a time tọ and at a later time t, respectively. The 
LMSD is related to the lateral diffusion coefficient, Dy, by 


D, = lim (AFijem(2)") 


t-00 4t 


(6) 


The pressure in confined systems by parallel plates is 
divided in the parallel and perpendicular directions. The 
parallel pressures, Pı, was obtained from 


P= 0.5(0xx F Oyy), 


where cx and oy, are the normal stresses in the x and y directions. 
The system structure was analyzed with the lateral radial 
distribution function (LRDF) gi(r), defined as** 


aO =a Xt- [O(la 51 +) - (lv -al -$)]. 


ify 


where 6(x) is the Dirac ô function and the Heaviside function 
6(x) restricts the sum of the particle pair in the same slab of 
thickness òz = ø. The lateral radial distribution function is 
proportional to the probability of finding a particle at a distance r 
from a reference particle inside the slab of thickness ôz. 

In order to check if the Janus system shows density anomaly 
we evaluate the temperature of maximum density (TMD). Using 
thermodynamical relations, the TMD can be characterized by 
the minimum of pressure versus temperature along isochores, 


The fluid, micellar and aggregated regions in the PT phase 
diagrams were defined analyzing the snapshots of the systems, 
the lateral diffusion coefficient, Dı, and the lateral radial 
distribution function, g(r). To define the nanoparticles in the 
same aggregate we defined a minimal distance equals fmin = 1.2 
and used a purely geometrical criterion: if the distance between 
one monomer of one dimer and a monomer of a distinct dimer 
is smaller than 7pin, then both dimers belong to the same 
cluster. 


Results and discussion 


We start our discussion showing the distinct micelles observed 
for the confined system. Regarding the temperature and density, 
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Fig. 2 Aggregates observed in our simulations: (i) trimeric cluster with 
Ne = 3, (ii) tetrahedral cluster with ne = 4, (iii) hexahedral cluster with n. = 6, 
(iv) spherical cluster with Nne = 8, (v) spherical cluster with Nne = 19, 
(vi) elongated cluster with ne = 10, (vii) elongated cluster with Ne = 20, 
(viii) disoriented rippled lamellae and (ix) oriented rippled lamellae. Blue 
particles are the A monomers and red the B monomers. All figures are 
projections in the xy plane. 


the Janus dimers aggregate in clusters with a distinct number of 
nanoparticles per cluster, ne. At lower densities, trimeric clusters, 
with ne = 3 nanoparticles in each aggregate and tetrahedral 
clusters, with n. = 4, are more common. Increasing the density, 
hexahedral clusters (n. = 6), as well as spherical and elongated 
micelles with distinct ną, are observed. The shape of each aggre- 
gate is shown in Fig. 2. For densities up to a threshold a 
coexistence of two or three of these micelles was observed. Above 
the threshold, one single rippled lamellar cluster with all the 
nanoparticles is observed. This lamellar phase can have a disor- 
iented structure, as shown in Fig. 2(viii), or an oriented structure, 
Fig. 2(ix). 

The region of micelle coexistence and the threshold to the 
lamellar phase depend on temperature, so we take the example 
of T= 0.10. In Fig. 3, we show the mean number of dimers in 
each cluster, (ne), as a function of the system density. We can 
establish a relation between (ne) and the type of aggregate. For 
0.05 < p < 0.10, (ne) varies between 4 and 5. Analyzing the 
system snapshots, for these densities we can see a coexistence 
of trimeric, tetrahedral and hexahedral clusters. As p increases, 


0.5 


Fig. 3 Mean number of dimers in each cluster, (ne), as a function of 
density for T = 0.10. Inset: Zoom in the region p < 0.30. Error bars smaller 
than the data point are not shown. 
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more hexahedral aggregates and less trimeric aggregates are 
observed in the solution. This region, with the coexistence of 
trimeric, tetrahedral and hexahedral clusters, was labeled as 
region I in Fig. 3. For the densities inside region II we observe a 
mixture of tetrahedral and hexahedral clusters, in region III 
hexahedral and small spherical clusters and in region IV a 
coexistence of small spherical and elongated clusters. All these 
previous aggregates where observed in the bulk simulation, 
with a similar (n,).°* However, the confinement disturbs the 
self-assembly as the density increases. Therefore, if for lower 
densities the system aggregates in the same micelles observed 
for the bulk case, for higher densities new kinds of self- 
assembled aggregates should be induced by the confinement. 
In this way, the spherical and elongated micelles can have a 
higher (ne) compared to the bulk case. Basically, the limited 
space induced by confinement leads two or more smaller 
micelles to merge into a large cluster. This region, with spherical 
and elongated micelles formed by more nanoparticles than observed 
in bulk, was labeled region V. The size of the micelles grows 
continually up to p = 0.35 for T = 0.10, with approximately 10 
aggregates with 100 nanoparticles each, as shown in Fig. 3, and 
above this threshold all the particles aggregate in a lamellar 
structure, region VI-A of Fig. 3. For temperature T = 0.10 and 
p > 0.40, the rippled lamellae are disoriented, without a 
preferred direction, as shown in Fig. 2(viii). However, for some 
values of temperatures and densities, as T = 0.25 and p = 0.375 
to p = 0.425, the lamellar phase has directional ordering, as 
shown in Fig. 2(ix). The region where we observed the oriented 
rippled lamellae was labeled VI-B. These lamellar structures, 
usually observed in amphiphilic molecules as in Janus 
particles,** °° were not observed in bulk for our model of Janus 
nanoparticles. Therefore, this new structure is induced by 
geometrical confinement. While in the bulk the particles do 
not have any geometrical restriction, remaining in the micellar 
phase when the density is up to p = 0.50,** the confinement, 
associated with a high density, leads the dumbbells to aggregate 
into a lamellar cluster. The time evolution of distinct lamellar 
phases is shown in Fig. 4. For lower temperatures, the entropic 
contribution to the free energy is not sufficient to change the 


t=0 t = 25000 


covet A ce} 
Li) 


Fig. 4 Lamellar phase for the density p = 0.40 at temperatures (A) T = 0.10 
and (B) T = 0.25 at three distinct times: end of equilibration time 
(t = 0), half of the production time (t = 25000) and end of simulation 
(t = 50 000). The figures are projections in the xy plane. 
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Fig. 5 Pı x T phase diagram for the confined Janus dimer system. The 
distinct micellar and lamellar regions discussed in the previous figures are 
indicated and separated by dashed black lines. The thin black lines are the 
isochores. The density anomalous regions is defined by the TMD line, and 
the diffusion anomaly region by the lines of maxima and minima diffusion. 


initial configuration and the lamellar structure does not change 
with time, remaining disoriented, as shown in Fig. 4(A). However, 
for higher temperatures, the initially disordered configuration 
changes to the oriented rippled structure, as we can see in 
Fig. 4(B). This lamellar phase with a preferable orientation is 
characteristic of dumbbell systems with one or two monomers 
that interact by a two length scale potential,?®?? but at lower 
temperatures it is disturbed by confinement. 

In Fig. 5 we show the qualitative P; x T phase diagram for 
the confined system. The distinct micellar and lamellar regions 
are indicated. As discussed previously, new self-assembled 
structures are induced by confinement. However, the aggrega- 
tion region does not shift to higher or lower temperatures. 
Curiously, at the temperatures we observed the oriented lamellar 
structure, for densities above p = 0.45, the system does not have a 
well-defined micellar structure, and an amorphous phase is 
observed. This phase is a fluid with small diffusion. Taking the 
isotherm T = 0.25 as a reference, the pressure drops when 
p > 0.45 and the system shows a reentrant fluid region. The 
plot of the parallel diffusion coefficient as a function of density 
for some values of T is shown in Fig. 6. As we can see that for 
T = 0.25 the system has Dı ~ 0 for p = 0.375, 0.40 and 0.425, and 
D, increases when p > 0.45, indicating a melting. The LRDF, 
shown in Fig. 7, also shows a decrease in the system structure, 
from the well-structured lamellar phase to the fluid phase. This 
melting induced by the increase of density was already observed in 
colloidal glass systems,”®”® but was not observed in our bulk system. 

The fluid phase also shows interesting properties, distinct 
from the bulk case. The first one is the density anomaly. As we 
can see, the isochores in Fig. 5 show a minimum. The dashed 
blue line, the so-called TMD line, connects these minimum 
points. The density anomaly was observed for pure anomalous 
monomers and dimers (AA dimers),”°** however was not present 
in the bulk Janus system.** The reentrant region, which occurs at 
higher densities, is located where the TMD line ends and splits the 
lamellar phase VI-B into two regions. 
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Fig. 6 Parallel diffusion Dı as a function of density p for different tem- 
peratures. The red lines indicate the minima and maxima in the diffusion 
coefficient. 
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Fig.7 Lateral radial distribution function (LRDF) g,(r) for T = 0.25, showing 
melting induced by the density increase. 


The bulk system has only diffusion anomaly, and it is 
preserved in the confined case. However, the diffusion extrema 
lines for the confined system, defined by the red dotted-dashed 
red line in the phase diagram, Fig. 5, have a temperature range 
larger than that in the bulk. Studies for confined anomalous 
fluids indicate that the TMD moves to lower temperatures for 
solvophobic confinement and higher temperatures for solvo- 
philic confinement.*°°?° For the bulk case, our hypothesis 
was that the TMD line was absorbed in the micellar region. 
Hence, since our confinement is solvophobic (the WCA potential 
is purely repulsive), it is surprising that the TMD appears, moving 
to higher temperatures. This leads to the question: why are 
anomaly lines shifted to higher temperatures? 

Gavazzoni and co-authors”? showed that anisotropy can 
shift the solid-fluid phase boundary of dumbbell systems made 
only by A monomers. More than this, the work argues that the 
kinetic energy, and, therefore, the temperature, has two contri- 
butions: translational temperature and non-translational 
temperature. Hence, the contribution from the dimer rotation 
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Fig. 8 Frontal (i) and side view (ii) of the system in the fluid region at 
T = 0.40 and p = 0.30. 


to the kinetic energy plays a significant role in this system 
behavior. In bulk, our Janus dumbbell can rotate freely in any 
direction, the only limitations are collisions with other dimers. 
However, confinement imposes a constriction to the system. 
In our strongly confined system, with L, = 4.0 and / = 0.8, not 
only translation in the z-direction is limited, but the combi- 
nation of confinement with the competing interactions of 
Janus systems lead to an interesting phenomenon. As shown 
in Fig. 6, the diffusion coefficient is distinct from zero — so, the 
particles are moving in the parallel direction. In Fig. 8(i) we 
show a frontal snapshot of the nanoparticle arrangement. As we 
can see, the particles are disordered - as expected for fluids. 
Notwithstanding, the side snapshot shows that in the z-direction 
the particles have a preferential position: the attractive B particles 
stay at the center, and the repulsive A particles are near the wall. 
Then, due to the confinement and Janus characteristics, the 
dimers are translated in the xy-plane, but without rotation. This 
places the A monomers side by side in the xy-plane, with an 
internal layer of attractive B monomers. This internal layer acts 
similarly to a solvophilic wall, with the B monomers pulling one 
another. As a consequence, the behavior is similar to water in 
hydrophilic confinement, and the anomalous region shifts to 
higher temperatures. 


Conclusion 


We report the study of confined Janus nanoparticles. This 
system has special interest in the design of new materials using 
confinement to control the self-assembled structures. We have 
found a rich variety of aggregates and micelles, including 
large micelles not observed in the bulk system. Moreover, two 
lamellar phases, with distinct orientations, were induced by 
confinement. The oriented lamellar phase region in the PT 
phase diagram is split by a reentrant fluid phase. This melting, 
induced by increasing the density, was not present in the bulk 
system as well. 

Another feature that was not observed in the bulk system is 
the density anomaly. The combination of confinement effects 
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with the competing interactions of Janus dimers shifts the TMD 
line to higher temperatures, exposing the anomaly that was 
hidden inside the aggregation region for the bulk case. More- 
over, the diffusion anomaly region increases, reaching higher 
temperatures. 

Our results show that materials composed by the association 
of a monomer which can be modeled by a two length scale 
potential and a standard LJ monomer show an interesting and 
peculiar behavior. New simulations are in progress to analyze 
the effect with distinct values of L,. 
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